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Abstract 

We study pion-nucleon scattering with a chiral lagrangian of pions, nucleons, 
and A-isobars. The scattering amplitude is evaluated to one-loop order, 
where Q is a generic small momentum, using a new approach which is equiv- 
alent to heavy baryon chiral perturbation theory. We obtain a good fit to the 
experimental phase shifts for pion center-of-mass kinetic energies up to 100 
MeV. A sigma term greater than 45 MeV is favored, but the value is not well 
determined. 
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I. INTRODUCTION 



The pion and the nucleoli play a central role in low-energy physics and there is a wealth 
of scattering data which allow us to test the application of effective Lagrangians. Several 
relativistic phenomenological models exist, which provide reasonably good fits to the 
experimental tcN phase shifts. In these models either the Bethe-Salpeter equation is solved 
approximately or the i^-matrix method is used to unitarize the tree amplitudes. Such 
models, however, do not offer a systematic approximation scheme. 

Chiral perturbation theory (ChPT) is a more attractive approach because it not 
only embodies chiral symmetry, which is fundamental to low-energy physics, but it also offers 
a systematic expansion in powers of the momentum. Further it ensures unitarity order by 
order. Gasser and Leutwyler [|] have shown that ChPT works nicely for mesons; however, the 
power counting fails when baryons are introduced 0. The power counting can be restored 
in heavy baryon chiral perturbation theory (HBChPT) ||10[ where the heavy components of 
the baryon fields are integrated out. Here we shall adopt a different approach which 
effectively removes the heavy fields after constructing the Feynman diagrams. This new 
approach preserves the power counting and gives results in agreement with HBChPT, at 
least to the order considered. 

Peccei []T2[ was the first to use a chiral lagrangian to calculate the vrA^ scattering lengths 
near threshold. Calculations with the more modern approach of HBChPT are discussed by 



Bernard et al. in a review and in Refs. ||Tj,|T5|. Particularly interesting is the recent 
calculation by Mojzis of the full amplitude to 0{Q^), with Q a generic small momentum 
scale. Mojzis employs just pion and nucleon fields. Since there are a number of unknown 
parameters, scattering lengths and effective ranges alone are not a stringent test of the 



approach. We mention that Datta and Pakvasa [jT^ have recently used the results of Mojzis 
to discuss low-energy scattering. However a fit to the phase shifts out to center-of-mass 
(cm.) energies in the A resonance region is desirable. This requires explicit consideration 
of the A field itself, along with the pion and nucleon fields. Here we shall carry out this 
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program by evaluating all diagrams up to one-loop order. 

The organization of this paper is as follows. In Sec. II we write down our effective 
lagrangian and establish our notation. In Sec. Ill we describe our method of separating 
out the soft contributions, which need to be calculated, from the hard contributions, which 
are subsumed in the parameters. Fermion loops and the implications for power counting 
are also discussed here. The formalism for ttN scattering is given in Sec. IV, with a listing 
of the expressions for the loop diagrams relegated to the Appendix. The calculated phase 
shifts are compared with the data in Sec. V, where we also discuss the effective ttA^'A'" and 
ttA^A couplings as well as the pion-nucleon sigma term. Our conclusions are given in Sec. 
VI. 



II. EFFECTIVE LAGRANGIAN 

Before writing down an effective chiral lagrangian we need to define our notation, much 
of which is fairly standard. The Goldstone pion fields 7r"(x), with a = 1, 2, and3 , form an 
isotriplet that can be written in terms of an SU (2) matrix: 

U{x) exp(2m(a;)/A) , (1) 

where /^r ~ 93MeV is the pion decay constant and the pion field is compactly written 
as 7r(a;) = tz{x) • |t, with being the Pauli matrices. The matrix U is the standard 
exponential representation and the "square root" representation in terms of ^ is particularly 
convenient for including heavy fields in the chiral lagrangian. The iso doublet nucleon field 
is represented by a column matrix 



N ^ 



( \ 

P 



(2) 



where p and n are the proton and neutron fields respectively. The A is a spin-| and isospin-| 
particle represented by an isoquadruplet field: 
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A. 



(3) 



A+ 

K 

It is convenient to introduce an isovector field A,, = TA,, in terms of the standard 2x4 



isospin I to ^ transition matrix: 



(4) 



where the isospin spherical unit vectors are cq = and e±i = Ti^x ± iey)/\/2. Explicitly, 
the components are 



A 



( v^A++^ 

A* 
AO 



(5) 



(6) 



(7) 



Following Callan et al. [0 we define a nonlinear realization of the chiral group SU (2)l ® 
S'f/(2)R such that, for arbitrary global matrices L G SU{2)i^ and R G S'f/(2)R, we have the 
mapping 



L®R: (e,iV,A^)^(e',iV',A;) 



(8) 



where 



i\x) = Li{x)h){x) = h{x)C{x)R^ 
N'{x) = h{x)N{x) , 



A' (x) = |/itr(/iT/iV)-A^(x) 



(9) 
(10) 

(11) 



As usual, the matrix U transforms as U'{x) = LU{x)K^. The second equahty in Eq. 
defines h{x) imphcitly as a function of L, R, and the local pion fields: h{x) = h{n{x), L, R). 
The pseudoscalar nature of the pion field implies h{x) G SU{2)y, with SU{2)y the unbroken 
vector subgroup of SU{2)i^ ® SU{2)^. The nucleon transforms linearly under SU{2)y as 
an isodoublet. While the isodoublet components of the isovector transform linearly 
in the same way as the nucleon field, the isovector itself is further rotated by the 0(3) 
transformation |tr(/iT/iV). 

Interaction terms invariant under the nonlinear chiral transformation may be conve- 
niently constructed in terms of an axial vector field a^{x) and a polar vector field Vf^{x) 
defined as 

i 11 

= + id.i^) =vi = \v,.T = (l - ^) 5.^] + • • • , (13) 

both of which contain one derivative. The polar vector field transforms inhomogeneously 
and the axial vector field transforms homogeneously: 



V 



hv^h) - ihd^h^ , (14) 



a'^ = /la^/i^ . (15) 

To maintain chiral invariance, instead of an ordinary derivative d'^, one uses a covariant 
derivative on the nucleon and A fields. These are defined by 

V^N = d^N + iv^N , (16) 
P^A^ = 9^ A^ + iv^r^u -v^x . (17) 

We also use the following definitions involving two and three derivatives on the pion field 

Vf,v = d^Vu - d^v^ + ^[v^,, v^] = -i[a^, a^] , (18) 
Df,au = d^ay + i[v^, a^,] , (19) 
DaV^y = d^v^y + i[va, v^y] , (20) 



all of which transform homogeneously in the same way as in Eq. ([T5|). 

To write a general effective lagrangian, we need an organizational scheme for the in- 
teraction terms. We organize the lagrangian in increasing powers of the fields and their 
derivatives. Specifically, as in Refs. P,p!9|, we assign to each interaction term a size of order 



g° with 

ft 

a = d+-, (21) 

where d is the number of derivatives on the pion field or pion mass {m.^) factors, and n 
is the number of fermion fields in the interaction term. That a is a characteristic of the 



interaction term is suggested by Weinberg's power counting ||20|, which we discuss in Sec. Ill 
below. Derivatives on the nucleon fields are not counted in d because they will generally 
be associated with the nucleon mass and not with the small momentum Q. In fact, as 
Krause has argued, it is 4^ — M that is of 0{Q^). Krause also counts 7^75 to be of 
0{Q^) and a single factor of 75 to be of 0{Q^). We adopt this counting for organizing the 
lagrangian, although we have argued that counting a single 75 factor to be of 0{Q^) is 
not precise. Our scheme allows a uniform organization of the pion self-interaction terms 
and those involving the heavy fermions. It differs from the "standard labelling" of Gasser, 
Sainio and Svarc where the number of fermion fields is not included. 

Taking into account chiral symmetry, Lorentz invariance, and parity conservation, we 
may write the lagrangian through quartic order [a < 4) as the sum of the order Q^, Q^, and 
Q'^ parts: 

C = C2 + Cs + C4 + AC , (22) 

where AC represents the counterterms, which can also be organized in powers of Q. We will 
adopt the counterterm method of renormalization, so we start with the physical masses and 
couplings and add the necessary counterterms. For simplicity, in this first investigation we 
will not explicitly identify the finite and divergent pieces of the various counterterms. The 
order part of the lagrangian is 
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£2 = iV(z^ + gA^l^a^ - M)N + \flii {d.U^d'^U) + \mlflii {U + - 2) 

+AlKbK + hA{A,-am + Na^.A,) + hAA;Yl^a,K , (23) 

where the isospin indices a, b = 1, 2, and 3, the trace is taken over the isospin matrices, and 
the kernel tensor in the A kinetic-energy term is 

A^'^' = -(if - M^)g^"' + i{-^>'V'' + YT^^) - l^i^Tp + Ma)Y , (24) 

suppressing isospin indices. Here we have chosen the standard parameter A = —1, because 
it can be modified by redefinition of the A field with no physical consequences In the 



ttA^A interaction of Eq. (|23| ) we have chosen the standard off-shell Z parameter to have 
the convenient value of — |. The value of Z has no physical significance since modifications 
can be absorbed in the other parameters in the lagrangian |23|. Similarly we simplify the 



ttAA interaction by choosing the offshell parameters defined in Ref. |2^ to be Z2 = — | and 
Zs = 0. 

With the notation 

5^= - (^^. -tv^,) , (25) 
we may write the order and Q"^ parts of C as follows: 

+ ^^iV7;. iVtr (a^a^) + gm^iViV tr ([/ + t/t - 2) + • • ■ , (26) 
£4 = ^mlN^,{U - U^)N + ^NYD^v,.N 

+^ml[K^.ti[id^{u - u^)t]n + mi[id^{u - u^)t]-a^} + ■■■ . (27) 

In Eqs. ( p6[ ) and (|27| ) the ellipsis represents terms that do not contribute to the vrA^ scattering 
amplitude. For example, the ellipsis in £4 includes the usual pion self-interaction terms with 
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four derivatives. The factors in £4 are introduced to correctly count the order of the 
symmetry breaking terms. However, it appears that counterterms of the above form with 
replaced by 5^ are needed even in the chiral limit [m.,^ = 0). Here 6 = Ma — M is the 
A-nucleon mass difference. We have applied naive dimensional analysis |p^ , p^ to the terms 



in Eqs. (|2B) and (|27|) so as to expose the dimensional factors. As a result, we expect the 



parameters to be of order unity. 



Using the pion and nucleon equations of motion [p0i^P6| , we have also simplified the 
contact terms listed in Ref. 0. For example, we reduce the 0{Q ) term N V^VuNtr [a^a^] 
to the sum of the 0{Q^) ki term, the 0{Q^) A4 term, and higher-order terms which we omit. 
As a result we have the minimum number of independent terms contributing to the vrA^ 
scattering amplitude up to 0{Q^). Note that the isoscalar-scalar and isovector- vector p 
fields given in Ref. |]19[ have been integrated out. Their effects show up in the contact terms 



Pn, and A2. For example, in terms of the pvrvr coupling {gpmr) and the pA^A^ coupling 
{gp), the rho gives a contribution to the A2 parameter of —2gpT,T,gpM'^f^/m'^. 

III. CHIRAL PERTURBATION THEORY WITH BARYONS 

A. Hard and Soft Contributions 

Given the effective lagrangian, one can derive the Feynman rules and carry out pertur- 
bative calculations of physical quantities in the standard way. However, as shown by Gasser, 
Sainio and Svarc the loop expansion no longer corresponds directly to the momentum 
expansion when we have heavy fermions. One way to overcome this difficulty is HBChPT 
||T0| , where the heavy components of the fermion fields are integrated out so that their effects 
on physical quantities only show up in the parameters of the lagrangian. Alternatively one 
can construct an explicitly nonrelativistic lagrangian pO] . 



We propose a different procedure here which involves manipulating the Feynman dia- 
grams themselves (see also the recent discussion of Gasser ||2^). First we obtain the Feynman 
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rules in the standard way. Then we separate the loop contributions into those from soft and 
hard momenta. We keep the soft contributions explicitly. These will have both real and 
imaginary parts in general, the latter being needed to maintain unitarity order by order. As 
for the hard contributions, we implicitly absorb them into the coefficients of the lagrangian. 
As we will see this procedure preserves a systematic power counting scheme. 

Specifically, we represent the hard momentum scale by the nucleon mass M. Other 
quantities of this order include the A mass Ma and the factor 47r/,r with 47r coming from a 



loop integral p8|. The soft momentum scale is denoted by Q, where Q <^ M. Quantities of 
this order are the pion mass m^r, the pion momentum, and the mass splitting 5 between the 
nucleon and the A-isobar. Also, as in HBChPT, we are interested in applications where the 
three-momenta of the external nucleons are of order Q. 

For the present we consider any loop diagram without fermion loops: we shall consider 
diagrams with fermion loops later. We obtain the unrenormalized soft part of the diagram 
by applying the following rules to the loop integral: 

1. Take the loop momenta of the pion lines to be of order Q. 

2. Make a covariant Q/M expansion of the integrand. 

3. Exchange the order of the integration and summation of the power series. 

Rule 1 ensures that the exchanged pions have soft momenta. As a result internal baryon lines 
will be nearly on shell throughout the diagram. Rule 2 implies that a covariant expansion 
of the baryon propagators in the integrand is made, which maintains the Lorentz invariance 
of the soft part. In Rule 3 the exchange of the order of integration and summation changes 
the result in general. Indeed the purpose of this maneuver is to remove the poles in the 
baryon propagators at hard loop momenta of order M. Clearly this is achieved because 
after application of Rule 2 the only poles in individual terms of the series are located in 
the soft momentum region of 0{Q). Of course the soft part obtained from our rules still 
contains ultraviolet divergences in the form of poles at c? = 4 in dimensional regularization. 
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\ 



FIG. 1. A one- loop nucleoli self-energy diagram. The dashed line represents a pion propagator 
and the solid line a nucleon propagator. 



We remove these divergences with the standard method of renormalization. Formally, if we 
denote the original integral by X, the unrenormalized soft part by SI and the renormalization 
operator by R, the final renormalized soft part is RSI. This is the physical loop contribution. 
As discussed in the next section, the soft loop contributions allow for a systematic power 
counting. 

As for the part of the original integral I that is discarded, namely I — RSI^ we call it 
the hard part. This hard part contains contributions from poles of the integrand at hard 
momenta. Thus, we should be able to absorb this part into the coefficients of the lagrangian. 



As Lepage |^9|] has argued from the uncertainty principle, large momenta correspond to 
short distances that are tiny compared with the wavelengths of the external particles, so the 
interactions must be local. 

At this stage a concrete example is useful. Thus, we evaluate the nucleon self-energy 
diagram shown in Fig. 0. From standard Feynman rules we obtain the self-energy: 



S.N(fc)-^./i J ^2^^. ' (28) 

where the free nucleon propagator 

and fi is the scale of dimensional regularization. We obtain the soft part of S^n by first 
making a covariant Q/M expansion of the integrand, while taking the pion momentum i to 
be soft. Since we assume k to be nearly on shell, we can make a covariant expansion of the 
nucleon propagator as 
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G{k + 



1 



2A; ■ £ + fc2 - M2 + 



26 



+ ■ ■ 



(30) 



2k-e + k^ - + te 

Here the contribution arising from first term in the square brackets is of order Q^^, while 
that due to the second and third terms is of order Q^^ x {Q/M). Subsequent terms will 
involve higher powers of Q/M. The precise way this expansion is carried out has to be 
tailored to the case at hand. For example, in vrA^ scattering we have k = {p + q) where 



and = . Then we may \ei k ^ p and £ ^ (£ + g) in Eq. (|30|) . 



Exchanging the order of the summation and integration in Eqs. (|2^) and 
the soft part: 



rM-3^?A. d'^^ 2k-i^-{}i + M)e 



+ ■ 



then yields 



(31) 



Here we explicitly show the leading contribution to the soft part S'S^^n which is of order 
Q^/M"^, as would be expected from Weinberg's power counting. The ellipsis represents 
higher-order terms. For illustrative purposes, it is useful to sum up the series of soft con- 
tributions. This can be carried out by noting that in the present case, after our exchange 
of the integration and summation, an i"^ in the numerator of any integrand can be replaced 
with in dimensional regularization. Then the exact soft part of the one-loop self-energy 
diagram is 



{2k-i + ml)^- {I/ + M)ml 



(32) 



Introducing 



u = 



2VA;2 



(33) 



which is of order Q, the integral of Eq. (p^) in the case < m,r can be written 



S^TT^{k) 



m 



2(47r^)nv^2 
A 



k^-M' 



A;2 



3gl le-M^ 



8(471^)2 



A;2 



(ml - 2u?)]i + 



2 

2u;m2 
~7W 



^ r~2 2 -1 

-- — \ mi — uo'^ cos 



-u 



+ M) 



327r2L + In ^ , (34) 



where 
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with Euler's constant 7 = 0.577- ■ ■. 

In the modified minimal subtraction (MS) renormahzation scheme, we obtain the renor- 
mahzed soft part ^SS^tn by including counterterm contributions (CTCs) to remove the term 
proportional to L in Eq. (Q) . We can further ensure that the pole of the nucleon propagator 
is at the physical mass with unit residue by additional on-shell mass and wavefunction coun- 
terterm subtractions as detailed in Eq. (|70D below. These CTCs can clearly be expanded 
in an infinite power series in (|(^ — M)/M. Thus, the divergences appear to all orders in 
the Q/M expansion. The divergences can be removed by introducing counterterms of the 
form j^^N{i<^ — M)"N with n an integer. However, we note that we may use the nucleon 



equation of motion p0|j25|j26[] to eliminate the above counterterms in favor of interaction 
terms involving multiple pions and nucleons. 

We may obtain the leading order Q^/M"^ contribution to Eq. (^) by approximating 
uj ~ (fc^ — M^) / (2M) and setting a/P ~ M in the denominators. It is not possible to expand 
this leading-order expression further since the square root and inverse cosine functions in the 
equation involve u and which are of the same order. Thus, we cannot absorb this soft 
part into the parameters of the lagrangian. This result is consistent with the expectation 
that the parameters should contain only high-energy contributions. 

The hard part of the self energy can also be evaluated directly. We find 



4 /2 ^ J {2tiY -M^ + ie){2k ■ i - - + ml + ie) 
-M^{M+ — -- — U ][ 2,271^ L + \a 



4(47r^)^' 



X 



32vr2L + ln^-l + V-^ . (36) 



Notice that the integral in the first equation of ( P^D is dominated by poles at momenta 
of 0{M). The final result can indeed be expanded in powers of {1^ — M)/M and thereby 

12 



removed by CTCs. If this were not done the power counting would be spoilt since the first 
term of the result is of 0{M) and the second is of 0{Q'^ / M), which can be contrasted with 
0(QVM2) for the soft part. 

While our procedure is plausible, a general proof would require consideration of diagrams 
of arbitrary complexity. Here we restrict ourselves to one-loop order for which it is easy to see 
that our procedure gives the same result as HBChPT. Indeed, the expansion of the baryon 
propagator in Eq. (^) generates the same effects as integrating out the heavy components 
of the baryon fields in HBChPT. For our example of the nucleon self-energy we can make 
the connection with HBChPT by introducing in Eq. ( ^I]) the four velocity such that 

= MVfj_ + with v^v'^ = 1. Projecting onto the light components by inserting the 
projection operators ^(1 = [^(1 fore and aft, and noting that 

1(1 +V) f75 +i) = 2i-S i(l +y) , (37) 

where 5^ = ^i750"^i/f we find 

2 



i(i +y) SE^k) i(i H) = 1(1 H) ^ V"'/ 



1 



+ ■■■. (38) 



V ■ q — V ■ I + te 

This is a well-known expression which follows directly from the Feynman rules of HBChPT 
|l^,[l6|. Our expression for the soft part, Eq. (|3^ , when projected onto the light components 



yields in leading order (or equivalently in the infinite nucleon mass limit) the result of 
Bernard et al. |T3|J3l obtained in HBChPT. 



Mojzis |16| has recently calculated diagrams for tiN scattering using HBChPT with just 
pions and nucleons up to one-loop order. The results of our calculation, given in the 
Appendix, agree with his (modulo differences in the parameterization of U and the treatment 
of finite terms arising from the product of a {d — A) factor with the l/(c/ — 4) pole, see below). 
For simplicity our examples have involved the nucleon propagator, but we emphasise that 
the delta propagator appearing in loop diagrams (denoted by an open box in Figs. ^ and ^ 
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FIG. 2. Sample diagrams with baryon loops not directly connected to external baryon lines. 



below) is treated in exactly the same way. This is necessary to preserve the power counting 



and it gives results which agree with the HBChPT approach of Hemmert et al. 32 



B. Fermion Loops and Power Counting 

Fermion loops were not discussed in the preceding section. We notice that there are no 
fermion loops in HBChPT after integrating out the heavy field components. Here we show 
that fermion loops have vanishing soft parts in our approach so that they can be ignored. 
First consider fermion loops that are not directly connected to external fermion lines, such 
as those in Fig. |^. We will work with nucleons for simplicity, although similar arguments 
can be given for A's. If we generalize our previous rules by taking the loop momentum to 
be of 0{Q), we can expand the propagators in the form 

«W = -Si(^+^'')(i + SJ + ---)- 

In dimensional regularization / d'^ii"- = so the soft part vanishes. Alternatively, a direct 
calculation of the loop integrals with the standard Feynman rules can be used to show that 
the contributions of the diagrams can be expanded in a power series in terms of small pion 
momenta {a. Q/M expansion) and can thus be absorbed in the lagrangian. In other words 
there is no soft part. 

Next we discuss fermion loops generated from four (or more) fermion vertices, such as 
those in Fig. |^. These can be considered to arise from Fig. ^ by shrinking the heavy boson 
lines to points. Since the momentum transfers through the boson lines must be small if it 
is feasible to integrate out the heavy bosons, the arguments given in connection with Fig. |^ 
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FIG. 3. Sample diagrams with baryon loops involving four baryon vertices. 




FIG. 4. Diagrams which yield Fig. 3 when the heavy boson propagators, represented by wiggly 
lines, are shrunk to a point. 

can also be applied here. Therefore the soft parts vanish. Alternatively, this may be shown 
directly by taking the loop momenta to be of 0{Q), using Eqs. (^) and (0), together with 
the relation 



where a is independant of i. Equation (^) is valid in dimensional regularization for any 



Finally, we need to consider the special case where an iV-baryon scattering process con- 
tains an A^-baryon intermediate state. A simple example is given in Fig. ^. We note that if 
we consider the four-fermion interaction in terms of the exchange of a heavy boson, as before, 
this diagram will not contain a baryon loop. As Weinberg noticed [^, there is an infrared 




(40) 



integers m and n (see Ref. |^0| for example). 



divergence in Fig. ^ for on-shell nucleons at zero kinetic energy. Indeed, the amplitude is 



proportional to 



P2W 
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FIG. 5. A two-baryon reducible diagram that needs special treatment. 



^2p2■^ + pl-M^ -le""'" ' ^^^^ 

where Pi{t) and P2{^) are polynomials in the loop momentum i. We have taken i to be of 
order Q and expanded the integrand in the manner of Eq. (|30[). The contour of integration 
is pinched between the two poles at £o = =t^e for Pi = P2 = (^,0), and so cannot be 
distorted to avoid these singularities. Of course, this just signals that our expansion fails. 
The way out of this difficulty has also been given by Weinberg: we should consider only 
iV-baryon irreducible diagrams for A^-baryon scattering processes. The reducible diagrams 
are summed up with the A^-body Bethe-Salpeter equation, the kernel of which is obtained 
from the irreducible diagrams. 

Even in one-baryon processes, such as ttN scattering, singular behavior can arise when an 
intermediate A goes on shell. A similar remedy is followed: first calculate the irreducible self- 
energy diagrams to a certain order, then sum up the string of reducible diagrams containing 
arbitary numbers of self-energy insertions (see Sec. IV). 

We can now discuss the power counting for irreducible diagrams that do not contain 
fermion loops. According to Eq. (|30| ) the leading order of a baryon propagator is 
and according to Rule 1 the loop momentum is of order Q. Thus all the power-counting 
arguments of Weinberg |^ carry over. It follows that the leading order of a Feynman 



diagram with L loops, Ej<( external baryon lines is Q" with 



16 



v = 2 + 2L-\E^ + Y.'^i(d, + \n, - 2) , (42) 

i 

where Vi is the number of vertices of type i characterized by baryon fields and di pion 
derivatives or m^r factors. (We used the quantity di + to characterize terms in the 
lagrangian in Sec. II.) In general each diagram may contribute at orders beyond the leading 
order. 



IV. PION NUCLEON SCATTERING 

We apply our formalism to vrA^ scattering and calculate the T matrix to 0{(^^\ The 
amplitude is obtained from tree diagrams constructed from our lagrangian £2 + -^3 + and 



one-loop diagrams constructed from £2- Following the standard notation of Hohler [33| and 



Ericson and Weise l^l] we write the T matrix as 

T,, = (TTfelTlTT.) = T+5,b + i[r,, rjT- , (43) 

where the isospin symmetric and antisymmetric amplitudes are 

T± = A± + i(^ + /)S± . (44) 

Here, as shown in Fig. ^(a), q and g' are the cm. momenta of the incoming and outgoing 
pions with isospin labels a and h respectively. The cm. momenta of the incoming and out- 
going nucleons are labelled p and respectively. The amplitudes and are functions 
of the Mandelstam invariants s = {p + q)"^, t = {q — q'Y-, and u = {p — q'Y- 



A. Tree-level Contact Terms & Nucleon Exchange 

In Fig. ^ we show the tree level Feynman diagrams arising from the contact terms and 
from one nucleon exchange, with the crossed diagram for the latter suppressed. The vertex in 
Fig. ^(a) arises from any of the interactions in £3 and £4 (except for the Ai and Xq terms), 
as well as the Weinberg term —N'-y^v^N . It is straightforward to obtain the amplitudes 
arising from the contact terms. The results are 
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(a) 



(b) 



FIG. 6. Tree-level diagrams for vrA^ scattering, (a) contact interactions; (b) nucleon exchange 
with the cross diagram suppressed. 



A; 



B? 



1 



1 



2(3^(2ml - t) - AK^ml + 
(ki - 2A4)(s - u) , 



8M2 



u 



[s -u] 



16M2 



u 



(45) 

(46) 
(47) 

(48) 



The parameters here will absorb the divergences arising from the one-loop diagrams. They 
depend on the scale of dimensional regularization, /i, in such a way that the complete T- 
matrix is //-independent. The contributions from nucleon exchange shown in Fig. ^(b) are 



well-known, see for example Ref. |3^. Including the crossed diagrams, we have 

4+ — ^ A/f 

J n 

9 A „^2 



J n 
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M2 



+ 



2/2 /2 VS-M2 M-M2 



(49) 

(50) 
(51) 
(52) 
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FIG. 7. One-loop A self-energy diagrams. The open box represents the free A propagator. 



Here we have used the exact nucleoli propagator, although it could be expanded in chiral 
orders as in Ref. |]T^. The difference would appear beyond 0{Q^) which is the level of 
precision of the present calculation. Note that identification of what is to be included to 
0{Q^) is model dependent since the hard momentum scale could be M, Ma, the average 
mass M, or 47r/^. Results obtained with different choices will differ at 0(<5^) and beyond, 
but this should not affect the quality of the fit to the data so that we shall simply write our 
expressions in convenient form. 



B. A Exchange 

When the A appears as an intermediate state for nN scattering, the tree-level T-matrix 
diverges at s = so the power counting fails. As argued earlier, we expect the power 
counting to work only for irreducible diagrams. Thus we evaluate the one-particle irreducible 
self-energy diagrams which, to one-loop order, are those of Fig. |^. Diagrams containing one 
or more self-energy insertions are then summed to replace the free propagator by the dressed 
propagator which is finite. 

We start by writing down the free A propagator 



GUk) = ^Zr^^^ [-9.^ + 17.7. + - + . (53) 



This may be recast in terms of the spin projection operators p5| , |36 



as 



where 
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1 2 

1/2 1 1 



1 

22 ~ 'i^'^^i'^v 



The spin projection operators obey the orthogonahty relations 

The dressed A propagator contains any number of irreducible self-energy insertions: 

IJ.U — + ^ fi\^ ^au + '^fj.X^ ^ar^ -kv ^ 5 l<3 ' J 

where S^j,(A;) is the A self-energy. Since ^^G^^^ik) and A;^G°^(A;) do not contain a pole at 
= M\, we conclude from Eq. ( p?] ) that the tensor terms in Tj^y{k) constructed from 7'^ 
and generate non-pole terms in the dressed propagator G^^,. It is not hard to see that 
these terms start to contribute to the T-matrix at order . We can thus greatly simplify 
our calculations by noting that in the A self-energy tensor 

= Sa + , (58) 

the aforementioned tensor terms, represented by the ellipsis, can be neglected. Renormal- 
izing Ea in the MS scheme we obtain S^^. We make additional on-shell mass and wave- 
function counterterm subtractions such that, when the imaginary part of the self-energy is 
neglected, the pole of the A propagator lies at the physical mass with unit residue. Thus, 
the final renormalized self-energy is 

where 3? refers to the real part. Breaking the A self-energy into real and imaginary parts 
we obtain the dressed A propagator 
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^ + Ma 



(p3/2) 



l/2^ 



21 )^J.l' 



mi 



2,, + Ma)(P2W- 



(60) 



We have ignored 0{Q^) contributions to the non-pole terms which involve P^/^ since, as 
we have remarked, they do not contribute to the vrA^ T-matrix until 0{Q^). Also in the 
numerator of the pole term in Eq. ( pUD we have neglected terms of 0{Q^), which contribute 
to the T^N T-matrix at 0{Q^). It is convenient to write the real part of the A "polarization" 
as a function of rj^ = — M'^)/M, with the mean baryon mass M = |(M + Ma). To 
0{Q^) the two diagrams of Fig. ^ yield 

4 hlMA 



3 (47rA)2 



25 h\MA 



Vl-ml)j{vk)-{s'-ml)j{6) 



iVk - S f - ml J{r]k - 6) 



m^(l + ln^)(%-5)} 



(61) 



where the function J is defined in the Appendix. For the energies considered here the 
imaginary part arises from the first diagram of Fig. ^ which gives 



TA{k' 



TT 



/)2 



x^(F-M2) 



k' + M' + 2MMa (k' - M'Y - {k' + 3M')mi 



2 — AkP'm'i 



(62) 



We have evaluated Fa up to 0{Q'^) because this significantly improves the accuracy of the 
A decay width, Fa(A;^ = M^). In fact the error is negligible when compared with an exact 
evaluation using tree-level coupling, g-^Ni^ = hAM/ (for the coupling to 0{Q'') see Eq. ( ^Sf) 
below). Specifically this is ||33| 



2 I 1 3 

■^exact _ 9ttNA \1 I 



where \q 



{4M^ -ml){S' -ml)/{AMl). 



M + J\q\'^ + 



(63) 
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FIG. 8. The A exchange diagram with the dressed propagator represented by a soHd box. 



Using the dressed A propagator of Eq. (^), the A exchange contribution to the real 
part of the T-matrix, pictured in Fig. is 



T, 



ba 



^ ^QMlAR{s)s[ai{s,t) + a2{s,t) 



(s-M^ + ml' 



M{s - M^) + (M + 2M^)ml + (s - + 2MMa + ml) 



(64) 



where we define 7] = \[s — M'^)/M and the real part of the A propagator is 



An{s) 



s-MI-Ua{v) 



(65) 



[s-Ml-U^{v)]^ + MiTl{s) 
We do not expand the propagator for the reason given at the end of Subsec. A. The 
definitions of the functions ai and a2 are 

(s- + ml' 



ai{s,t) = 2M 



ml - it - |(s-M' 



6s 



M{s - M^) + ml{M + 2Ma) 



a2(s,t) = - \t+%MM 



[s - + ml] 
6s 



(66) 



To this should be added the result for the cross diagram which is obtained from the above 
by the replacement s u and the interchanges a ^ b and q — g'. 



C. One-Loop Diagrams 

A set of one-loop diagrams that contribute to the vrA^ T-matrix at 0{Q^) is shown in 
Fig. ^. Here we can use the free A propagator since no singularities are generated in the 
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FIG. 9. A set of one-loop diagrams which contribute at 0{Q^). Crossed diagrams for (d) to 
(n) are not shown. 
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(a) (b) 



(c) 



(d) (e) 



(f) 




FIG. 10. Diagrams with one-loop vertices which contribute at 0{Q'^). Crossed diagrams are 
not shown. The sohd circle in (a) and (g) refers to Ai and Xq vertices, respectively. Each diagram 
implicitly includes its counterpart where the lower vertex is dressed. 



T-matrix. A covariant expansion of the A propagators, as well as the nucleon propagators, 
is made in the manner discussed in Sec. IIIA so that the denominators are of order Q. These 
A propagators are denoted by open boxes in Fig. ||. 

As an illustration, we calculate the T-matrix for Fig. ^(e). Using the standard Feynman 
rules, we find 



16/4 ^ J {27rY e^-ml + %e 

xu{p')l^,T'G{p' + m'i,T'>G{p + q + mi^T'^Gip + £)f75r%(p) . (67) 

Noting that p'^ = Pfi + ^ I'fj,, and that q, q', and £ are of order Q, we can expand the 
nucleon propagators as in Eq. ( pO| ) and select the leading terms. Then the contribution to 
0{Q') is 



d 



16/4 - J (27r) 

ujp'W - M)iy + M)i(^ - M)f((5"'' - |[r^ r-])uip) 
^ (£2 -ml + ie){2p ■ £ + iey{2p ■ i + 2p ■ q + ie) ^ 
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As in Sec. Ill, we can cast this in HBChPT form. Setting = Mv^ for the incoming 
nucleon, we easily obtain 



u{p'){l ■ Sfjq' ■S){q- S){6^' - |[r^ T^])u{p) 
^ {& -ml + ie){ve + ie)^{ve + v-q + ie) ' ^ ^ 



It is tedious, but straightforward to evaluate the numerators of Eqs. (^) or (|69D and carry out 
the integration using the results in the Appendix. We use the modified minimal subtraction 
scheme (MS) to remove the singularities and the resulting real part of the T-matrix is listed 
in Eq. ([A20|) . Note that divergent terms proportional to l/{d — A) involves polynomials in 
the variables and so they can be absorbed into the contact-term contributions. For the same 
reason, we can absorb finite polynomial terms obtained from the product of a l/{d — 4) 
singularity and {d — 4) factors. The latter may be generated from d- dimensional 7-matrix 
algebra or from integration of tensor factors such as i^i^, in the numerator of the integrand. 
The above simplification has been exercised for the results listed in the Appendix. 

The real parts of the T-matrix for the remaining diagrams of Fig. ^ are evaluated in like 
fashion to 0{Q^) and listed in the Appendix. Note that the diagrams in Fig. ^(m) and (n) 
involve the nucleon one-loop self-energy which is real for the energies of interest here. As 
with the A, we renormalize in the MS scheme to obtain and make additional on-shell 
mass and wavefunction counterterm subtractions. Thus, the renormalized self-energy is 

ETik) = S^(fc) - E^ik) 1^^^^ 1^=^, i^-M). (70) 

Since these diagrams do not give singular contributions to the T-matrix we do not sum the 
self-energy insertions. 



We also need to evaluate the vertex modification diagrams shown in Fig. |10[ The A 
propagators denoted by open boxes have the same meaning as for Fig. |^. However there are 
now cases where the A can go on shell for which we use the dressed propagator of Eq. (^0|) , 



as indicated by the solid box. For the diagrams in Fig. |TO|(i) and (j) the contributions 
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FIG. 11. One-loop diagrams which either vanish, (a) to (f), or do not contribute at 0{Q^), {g 
to (1). Crossed diagrams and diagrams with the time ordering reversed are not exphcitly shown. 



proportional to Ar(s) vanish at the pole where s = M\ (and so ?7 = 5). As a result, at the 
pole there are important contributions to the T-matrix coming from the imaginary part of 
the nNA vertex. These contributions actually are of 0{Q^) for energies very close to the A 
resonance {\ri — 5| ^ w^tt)- Thus, we keep the leading-order pole corrections in Eqs. ( |A41| ) 



and ( |A42| ) . We have also considered the next-to-leading order {0{Q'^)) pole corrections for 
Fig. |TD|(i) and (j), as well as the leading 0{Q^) pole correction to the A exchange diagram of 
Fig. 1^. These gave a very small effect and were not included in the final fit. The expressions 
for the vertex modification contributions to the T-matrix are also listed in the Appendix. We 
reiterate that our expressions for diagrams that do not involve the A propagator agree with 



those recently given by Mojzis |16]. Finally we show in Fig. O a set of one-loop diagrams 



which do not need to be considered. The top row, |TT](a)-(f), is identically zero due to the 
vertex structure or to violation of isospin conservation. The lower row, [Tl](g)-(1), does not 
contribute to 0{Q^). This is a welcome simplification. 

From the real part of the T-matrix the real part of the elastic scattering amplitude, /„, 



is obtained by the standard partial wave expansion Here the isospin-spin partial wave 
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channels are labelled by a = (/, 2/, 2 J) with / the orbital angular momentum, / the total 
isospin, and J = / ± | the total angular momentum. The phase shifts 5a are then given by 

^fa = 1 — r sin (5q, cos (5q, . (71) 



D. Pion— Nucleon a Term 



We may obtain the nucleon a term from the Feynman-Hellman theorem, 

cr(0) = ml 



Working out the nucleon self-energy up to one-loop order, we find 

° M " 2(47r/,)2 " 



(72) 



+ 



8 h\ 



3 (47r/0 



{5'-mi)J{5) + \mi5\Yi^ 



(73) 



where the function J is defined in the Appendix and Mq is the "bare" nucleon mass defined 
to be independent of the pion mass. Note that, as in HBChPT, contact terms proportional 
to 5^ are needed to absorb the divergences and /z dependence associated with the A degrees 
of freedom. We have not specified these explicitly, hence our Mq depends on 5 and /x. The 
a term up to 0{Q^) is then 



a(0) 



4k2 2 9 vr^i 3 8 h\ml 
-ml — -— — :r^m^ ' 



(74) 



4(47rA)2 - 3(47r^)2 

Note that the /i dependence in J{6) may be absorbed in K2 so that a is independent of fi. 

The isospin even on-shell amplitude D^{iy,t) at the Cheng-Dashen point can be related 
to the irN a term, a{2mD, as follows ||38[| : 



f^D^iu = 0, t = 2mi) = a{2mi) + 0{ml) 



where v = {s — u)/ (4M) and D is obtained from 



(75) 



L)+(z/, t) = A+{u, t) + z/5+(z/, t) 



(76) 
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by the subtraction of the nucleon pole term which comes from the tree diagram in Fig. ^(b) 
and the one- loop diagrams in Figs. |10|(a) to [10|(f). At the Cheng-Dashen point only Figs. 
^(a), |^(b) and Figs. ||(r) to ^(u) contribute to the a term. Using the results for the T-matrix 
as given in the Appendix, we find 



a{2mi) - a{0) = - 



3 rrgl 



3 , KK 



4 (47r^ 



3MIS 



2 hlml 

I A. TV 



3 (47rA)2 
+ / dx- 



(tt- 4)5 -4^52 



ml In I 



5^ 



252 



:ln ' • V • ' . (77) 

/52 - m2(l - 2x + 2x2) 5 - ^b'^ - ml{\ - 2a; + 2x2) _ 

Note that since we count b to be of order Q the second term on the right is of the same 

order as the other terms. Equation (|77D agrees with the result of Bernard et al. [p^^ who 

obtain this expression from HBChPT and dispersion relations. (Their integral is in a form 

which is slightly different from ours, but both integrals give the same numerical result.) 



E. ttNN and vriVA Couplings 

The nNN vertex up to one-loop order consists of the tree vertex generated from the axial 
term in the lagrangian and the one-loop diagrams shown in the upper part of Figs. p!o|(a) 
to (f). Following our procedure, we can straightforwardly calculate the one-loop vertex 
function r°(fc, k', q), where k {k' ) is the incoming (outgoing) momentum of the nucleon and 
q = k' — k is the momentum transfer. The irNN coupling for on-shell nucleons is obtained 
by sandwiching this vertex function between nucleon spinors: 

u{k')T%k, k', q)u{k) = g.NN{q')u{k')j,T^u{k) . (78) 

We find 

(0) + 0{q'Q') , (79) 

which implies that the difference between the on-shell pion coupling g-^NN^m^) and gnNN{0) 
is of 0{Q^). This justifies the usual assumption of a small variation in gnNNiQ"^) between 
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= and = when the Goldberger-Treiman relation is derived To 0{Q'^) we 



obtain 

gnNN{0) = < 1 



mz 



gA M2 3 (471/^)2 - /i2 ' 6 (AnUy 
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In ^ + 
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275 (47rA)2 
200 h\hA 



irmi + {6' -mi)J{6) + \mib In 



243^?A (47rA)^ 



25^-^2 -35J(5)-fm^ In ^ 



(80) 



Notice that Ai absorbs the divergences and //-dependence arising from the one-loop vertices. 
With the parameters obtained from fits to the ttA^ phase shifts Eq. ( |80D allows a test of the 
Goldberger-Treiman relation. 

Similarly, we can calculate the ttNA vertex from the tree-level term and the one- 



loop diagrams shown in the upper part of Figs. |TO|(g) to (1). Here, in the vertex function 
r^"(A;, A;', g) the label k now refers to the incoming A momentum. The nNA coupling is 
obtained by sandwiching this vertex function between the nucleon spinor and the A spinor. 



(81) 



We find 



(82) 



where 



gnNA{0) = — < 1 
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(83) 



35 (47r^)2 

Here \q plays a similar role to Ai for g^NN- The value of g^NA is complex because the 
intermediate pion and nucleon states for Fig. |l3(i) and (j) can go on shell. 
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V. RESULTS 



As fixed input parameters, we use the standard baryon and pion masses: M = 939 MeV, 
Ma = 1232 MeV, and = 139 MeV. We also take g = 92.4 MeV from cliarged pion 
decay, gA = 1-26 from neutron /? decay, and Ha = 1.46 from Eq. (|63D with gnNA = hAM/ f^^ 
and the central value of the A width = 120 ± 5 MeV. The physical results should be 
independent of the scale of dimensional regularization. To confirm that this is indeed the 
case we carry out calculations with three values of the scale, namely fi = Ma, = 1 GeV, 
and fi = 0.75 GeV. We then have eleven free parameters left: /J^r, /^tt! ^i, ^2, ^a, and Ai to 
\q. These are obtained by optimizing the fit of our calculated S- and P-wave phase shifts 
to the vrA^ scattering data of Arndt |^ . Because negligible error bars are given in the data 
at low energies, we assign all the data points the same relative weight in a least-squares fit. 
We fit the data for pion cm. kinetic energies between 10 and 100 MeV. 

The fit obtained for the S- and P-wave phase shifts as a function of the pion cm. kinetic 
energy is indicated by the solid line in Fig. |12[ Here the renormalization scale is chosen to 
be 1 GeV and the corresponding parameters are Set Al in Table 1. The experimental data 
points of Arndt |^ to which we fit are displayed as triangles in Fig. |12|; we also display 
there the older data of Bugg |^ (squares) and Koch and Pietarinen |^ (circles). The 
solid line shows a good fit up to 100 MeV pion cm. kinetic energy, slightly below the 
A resonance at 127 MeV. It is remarkable that a good fit is achieved, even with eleven 
parameters, since the T-matrix contains a number of complicated non-polynomial functions 
of the invariant variables. The agreement for most of the phase shifts extends beyond 
100 MeV, but 5531 starts to deviate markedly. Extending the range of cm. energies used 
for the fitting does not change the situation significantly. However, 0{Q^) contributions 
may become significant above approximately 100 MeV. This assertion is supported by test 
calculations in which we have included a few selected contributions of this order. However, 
without considering all 0(Q'^) contributions no definitive statement can be made. A better 
strategy is to examine the renormalization scale dependence of the phase shifts. The dashed 
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FIG. 12. The S- and P-wave phase shifts as a function of the pion cm. kinetic energy. 

The sohd, dashed, and dotted curves are calculated with parameters sets Al = 1 GeV), A2 

{fx = 1.232 GeV), and A3 {fx = 0.75 GeV) respectively. The data are from Arndt |^ (triangles), 
Bugg (squares), and Koch and Pietarinen (circles). 
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FIG. 13. The dependence of i sin 2^33 upon the pion cm. kinetic energy. The phase-shift data 



from Arndt |42] are indicated by the triangles. The sohd, dashed, and dotted curves are calculated 
with parameters sets Al, Bl, and CI respectively. 

and dotted curves in Fig. |12| display results obtained with /i = Ma and fi = 0.75 GeV 
respectively; the corresponding parameters are labelled Sets A2 and A3 in Table 1. The 
values for the three different values of the scale are similar and the phase shifts agree 
up to 100 MeV. Beyond this energy a dependence on the scale starts to appear indicating 
that 0{Q^) contributions are needed. The deviations are most noticeable for 6331 and 6pi3 
which suggests that they are the most sensitive to higher-order contributions. We remark 
that similar fits can be obtained with values of the renormalization scale as small as m^r. In 
Fig. [1^ there is a small gap in the 5p33 phase shift just below 100 MeV. Here unitarity is 
slightly violated and a phase shift cannot be determined. The point is made in Fig. [T^ where 
|sin25p33 is plotted for Set Al (solid curve). The maximum magnitude of | is exceeded in 
a 10 MeV energy region by < 10 %. This is a small violation and the salient point is that 
sensible results are obtained without resorting to a phenomenological i^T-matrix approach to 
enforce unitarity. 

While we expect the parameters of Table 1 to be natural, i.e. of order unity, some of 
them are close to 10. There could be several reasons for this. First, we have no definite 
prescription for assigning factors of 2 in the definition of the parameters. Second, the phase 
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shifts are derived from the T-matrix where there are strong cancellations among the various 
terms. This suggests that the parameters have significant uncertainties. Third, we note that 
the parameters are scale dependent. They appear more natural for Set A3 with fi = 0.75 
GeV than for Sets Al and A2. An example of this is the parameter Ag which contributes to 
the effective vrA^A coupling, see Eq. (|55D. In fact the last row of the table shows that the 
effective vrA^A couplings obtained from sets Al, A2, and A3 are quite similar and show a 
very small deviation from the tree-level value. Thus, we should not be unduly concerned by 
the size of some of the parameters. 

In Table 1 we also show the value of the effective tcNN coupling, which remains close to 
the tree-level value. The ttNN coupling compares quite well with the determination of Arndt 
et al. ^3 of Q-knn / {^^^) = 1-03. We also tabulate the nucleon sigma term calculated from 
Eq. (0). Essentially the same value is obtained for Sets Al, A2, and A3 so the sigma term 
is renormalization scale independent, as it should be. The magnitude of 108 MeV, however. 



is much larger than the generally accepted value of 45 ± 8 due to Gasser et al. . It is 
only slightly reduced to 92 MeV when we fit to the older data of Koch and Pietarinen [HB . 
The value of cT(2m^) — (t(0) = 15.7 MeV depends only on qa and and is uncontroversial 
I TB| , p!5| , ^ . Thus our predicted (j(2m^) is also much larger than the value of 64 ± 8 MeV 



obtained using dispersion relations. In order to see how sensitive our fit is to the sigma 
term, we make two other fits in which (t(0) is constrained to be 75 and 45 MeV respectively. 
This is achieved by fixing ^2 via Eq. ([7^), while allowing the remaining parameters to vary. 
We choose the scale of dimensional regularization to be /i = 1 GeV. The parameters thus 
obtained are respectively labelled Bl and CI in Table 1. The corresponding phase shifts are 



denoted by dashed and dotted lines in Fig. 14 and these can be compared with the solid 



curve obtained with Set Al. The value of for Set Bl is 50% larger than for Set Al and 
only the ^531 phase shift shows any visible difference between the two cases for cm. energies 
below 100 MeV. The tcNN and ttA^A couplings are also similar. However, for set CI 
is more than a factor of 2 greater than for Set Al and larger differences appear for the 
phase shifts, particularly the S"- waves, and couplings. Clearly a{0) is not well determined, 
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FIG. 14. The S- and P-wave phase shifts as a function of the pion cm. kinetic energy. The 
sohd, dashed, and dotted curves are calculated with parameters sets Al, Bl, and CI respectively. 
The data are as in Fig. 12. 
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although values larger than 45 MeV are suggested. Note that Mojzis [0 obtained a value 
of 59 MeV by fitting to the threshold parameters, the ttNN coupling, and the sigma term. 
The parameters for Sets Al, Bl, and CI show some significant differences indicating, as 
mentioned earlier, that they are not well determined by fitting the phase shift data alone. 
We also note that Fig. |T^ shows that Sets Bl and CI violate unitarity at about 150 MeV 
cm. energy, although here, as we have remarked, we expect 0{Q'^) effects to be important. 
However, it may be significant that just below 100 MeV Set CI gives a larger violation of 
unitarity than the other parameter sets. 

It is also interesting to examine the threshold (vanishing pion kinetic energy) S-wave 
scattering lengths (ci2/) and the P-wave scattering volumes (ci2/2j)- The experimental values 
were not used in the fit and they are compared with our predictions with Sets Al, Bl, and 
CI in Table ||. The agreement for the P-waves is reasonable with not too much sensitivity 
to the parameter set, although 033 does become rather small for Set CI. For the S'-waves it 
is instructive to examine the isoscalar and isovector scattering lengths, Bq = (ai + 2a2)/3 and 
bl = (0.3 — ai)/3. For 61, Ref. |^|, Ref. [Q, and a recent determination by Sigg et al. 
give —0.091, —0.091 ± 0.002 and —0.096 ± 0.007 respectively. These are in agreement and 
somewhat larger in magnitude than our predicted value of —0.082, virtually independent 
of parameter set. The value of 60 from Ref. is consistent with zero, while Refs. |]46|j48[| 
give -0.010 ± 0.003 and -0.008 ± 0.007. Our parameter sets Al, Bl, and CI give 0.003, 
—0.013, and —0.023 respectively. From this result, and the others that we have discussed, 
we conclude that parameter set CI is less favored than the other sets. Note that we do not 
show the scattering lengths and volumes for Sets A2 and A3 because they differ negligibly 
from Set Al, indicating //-independence. We make the obvious remark that the accuracy of 
the predicted scattering lengths and volumes can be improved by including them in the fit 
at the expense of a poorer fit to the phase shifts. In the case of Set A this reduces the sigma 
term by about 10%. 
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VI. CONCLUSIONS 



We have introduced a new approach to chiral perturbation theory with baryons which 
involves manipulating the Feynman diagrams directly, rather than integrating out the heavy 
components of the baryon fields at the lagrangian level as in HBChPT. Our approach pre- 
serves Weinberg's power counting so that a systematic expansion in powers of a generic 
soft momentum scale, Q, remains valid. We can maintain relativistic covariance and apply 
the standard Feynman rules and 7-matrix algebra. At one-loop order we obtain equivalent 
results to HBChPT; it would be interesting to compare the approaches in higher order. 

We have calculated the T-matrix for nN scattering up to one-loop order. There were 
several components to this calculation. First, the chiral lagrangian introduced a number of 
contact terms with unknown coefficients which contributed to the T-matrix. Next, there 
were the well-known nucleon and A exchange diagrams. The former was straightforward, 
but the latter was singular at tree level when the A went on shell. The problem was solved 
by summing up 0{Q^) self-energy insertions in the propagator to all orders. The resulting 
amplitude then (approximately) obeyed unitarity for pion cm. energies up to at least 140 
MeV. Finally, we needed to evaluate the large number of one-loop diagrams which could be 
constructed with the tt, N and A vertices. 

We have performed a least-squares fit to the S- and P-wave phase shift data to determine 
the parameters involved. We were able to obtain a good fit up to 100 MeV pion cm. 
kinetic energy, slightly below the A resonance. Our predictions for the zero-energy scattering 
lengths were also reasonable. These results were independent of the renormalization scale, 

as they should be. Since r]/M is an expansion parameter and at the A resonance it 
becomes 5/M — 0.27, one might expect 0{Q^) effects to become important beyond the A 
resonance energy. This was supported by the scale dependence of our results in that energy 
region. Though straightforward in principle, it would be very tedious to include the 0{Q'^) 
corrections. The fit that we have obtained here contained an unavoidably large number of 
parameters, but it must be borne in mind that the T-matrix is a much more complicated 
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function of the invariant variables when loops are present than at tree level. Thus we 
regard the fit that we have achieved as a vindication of the basic approach. It seems, 
however, that one cannot determine the parameters with great accuracy using the phase 
shift data alone. In particular we were able to obtain very similar fits with parameters that 
corresponded to a nucleon a term of 105 and 75 MeV. We judge that the fit was slightly 
degraded when the currently accepted value of (j(0) = 45 MeV was adopted, so that a 
somewhat larger value appears to be favored from the analysis here. However this might 
indicate that O(Q^) effects need to be calculated in order to achieve an accurate one-loop 



result, as Ecker has recently argued |49]. It would be interesting to attempt to constrain 



the parameters further and to compare with the parameters found at zero energy ||T5|,|T6 



The latter would involve disentangling the A contributions as well as other effects that, for 
simplicity, we have absorbed in the contact terms. 

We thank S. Jeon and J. Kapusta for useful discussions and critical readings of the 
manuscript. We acknowledge support from the Department of Energy under grant No. 
DE-FG02-87ER40328. 



APPENDIX: 



We first list the integrals we need. Using dimensional regularization we obtain 

r ({d^ 2 1 

where L, defined in Eq. (pSj), is singular and the finite part 



m2 / Am? a/1 - ^ + 1 
lix) = 1 - In ^ - W 1 - ^ In for x < . (A2) 



The next integral to consider is 



{27fY {P - ml + ie){2p ■ i + z + ie) {Aixf^/p^ \2Vp^ , 

i 



IGvrp^ 



e{z)e{z' - 4p'mi)Jz^ - Ap^ml , (A3) 
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where 9 denotes the Heaviside step function and 



J(x] 



X — xln^ — \ — ml In ■ ^ 

7^ V ^ ^ y/^ 



(A4) 



X - xln^ - 2^ml - x"^ cos"^ ( - ^) (|a;| < m^r) 
Note that J(0) = — vrmTr. The next integral we need is finite: 

d'^i 1 1 



in 



4-d 



(A5) 



where 



/ Ami 1 
1 ^ sin"^ 

X 



4mJ. 



for x < 



(A6) 



Finally we require the finite integral 



{2ttY (£2 -ml + ie)[{i + Qf - ml + it]{2p ■ i + - Ml + ie) 



1 



327r2M5 



dx K{x) 



(A7) 



where 



K{x) 



In 



5+Y^(52-m2+x(l-a;) 



^ tan' 



((52-m2 +x(l-x)g2 > 0) 
{6^ - ml + x{l - x)Q'^ <0) 



(A8) 



i/m2-x(l-x)Q2-52 

In Ref. [|16| the expressions for the integrals of Eqs. ( |A1|) , ( |A3|) , and ( [X5| ) are written —Io{Q'^) 
— Jo{z/^/4j?)/y/4j^, and —Kq{0, Q^) / they agree with our results above. Our expres 
sions are written in terms of M, S, and rj and we recall the definitions for reference 



M = \{M + Ma) , (5 = Ma - M , 77 
It is useful to define the functions 



2M 



(A9) 



+ (r/-5)(252 



m^ln^ 



(AlO) 
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1 



riS{ri - 5) 
1 



(?7^ - ml)5J{r]) - (5^ - ml)r]J{S) + 7rml{r] - 5) 



m: 



1] — b 

These functions are finite, for example 



(All) 
(A12) 



77? 



Then the contributions from Figs. ^ and |10| to the real part of the T-matrix 
scattering are: 



9a 



9a 



96 



96 



9e 



rpba 
^9/ 



9 X 3(m2 - 2t)(2ml-t) ^, / 



12/2 (4^J^)2 

1 r 777.'^ 

—-.{ - v{8ml - 5t)I{t) + I0r]ml In -f 



9l 



Ami — t 

2 

7 



— Gtt 



2/7^ 1 



- M77 + i(27772 - t)J [2777^ + ^^47772 - t S{t) 

9/2 (47rA)2 L^^'^^' ~ ~ ^'^'^^^^^ + " ^^^^^^^ 

-(m^-2t)(252- 2777^ + t) / da;ir(x) - 1777^5 In ^ 

JO /7"' 

+ 



Mi - Mr] + i(277i2 - t) + 4775J J(5) 
(m^ -Mr] + \{2ml - t))^ + {25' - \ml + |t)7/]/(t) 
[Mi - M77 + ^(2777^ - t)) (52 _ ^2 ^ It) 



1 - 777.^ 



+27^5(52 _ ^2 ^ i^)j ^r^x[x) + fm^r/ln-f 



12/2 (4^J^)2 



(4777^ - t)J(t) - 2777^ In 



777: 



2 n 



fl (4vr/.)2 L 



r/J(r/) + §777^ In^ 



4/.^ (47rA)^ 



9/.2 (4vrA)^ 



X 



M^-M77-|(27772 -t) + 3r/^ 



(A13) 

for vrA^ 

(AM) 



(A15) 



(A16) 



(A17) 
(A18) 
(A19) 
(A20) 



(A21) 
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40 



rpba 
J-lQf 



rpba 
^109 



ba 
lOh 



pba 
lOi 



ba 
lOfc 



ba 
101 



T ,r 



,ab ,i[r^r"; 



5 



243/2 {A-kUYt^ 



Mr] + W)F,{r]-6, -6) 



3/2 M2 

In ^ 

9/2 (47rA)2 - /,2 



ai{s,t) + 0:2(5, t)^ 



r72 7,2 



9/.2 (47rA)2 



MaTUs) 



8h\ 



1 



81/2 (4^J^)2 



a;i(s,t) + 0^2(5, t)^ 



8W (47r/.) 
200 /i2, /i2. 



ai(s,t) + 02(5, t)^ ^3(77 - 6, -6)Ar{s) 



(A38) 
(A39) 
(A40) 



(A41) 



(A42) 
(A43) 
(A44) 



729 /2 (47rA)2 

Notice Eqs. (A33) - (A36) for diagrams 10(g)-10(l) have been written in terms of ai, 02 of 
Eq. ( 166|) rather than truncating these expressions to OIQ"^). This is necessary to ensure that 
the resonance only contributes to the P33 channel and does not "contaminate" the other 
channels. To the above contributions should be added the results for the cross diagrams 
where applicable. They are obtained from the listed expressions by the replacement s — > m, 
7] ^ fj = {u — M2) / (2M), and the interchanges a ^ b and q ^ — g'. 
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TABLE I. Parameter sets obtained from fitting the ttN phase shifts for various values of the 
renormahzation scale /x. The deduced TrATiV and TriVA couplings and the nucleon a term are also 
given. 





Al 


A2 


A3 


Bl 


CI 


M(GeV) 


1.00 


1.232 


0.75 


1.00 


1.00 




-3.5790 


-4.2084 


-1.9714 


-3.1385 


-3.1641 




2.3219 


3.5370 


0.1993 


1.9697 


2.3169 


Ki 


5.6578 


6.5464 


1.5485 


4.6914 


5.2765 


K2 


-2.3683 


-2.5845 


-2.1417 


-1.993lt 


-1.6286t 


Al 


-2.6652 


-1.9338 


-4.2176 


-2.2593 


-9.7605 


A2 


-6.4822 


-4.5533 


-7.7712 


-8.2178 


-11.7498 


A3 


11.0504 


11.9228 


4.6877 


10.7125 


15.2134 


A4 


-4.4304 


-2.5495 


-6.3775 


-6.1955 


-9.6040 


A5 


4.9702 


3.5669 


4.1456 


5.3523 


8.0874 


Ae 


9.3759 


9.2927 


4.8363 


7.3192 


7.4099 


h-A 


1.6243 


1.4418 


1.0986 


1.6779 


2.4626 


a(0)(MeV) 


106 


109 


108 


75 


45 




1.0445 


1.0378 


1.0409 


1.0091 


0.9585 


\9.NA\/C-i^) 


0.9898 


1.0134 


1.0337 


1.0532 


1.1294 


^ Fixed so as 


to obtain the listed value of cr(0). 









45 



TABLE II. The calculated S-wave scattering lengths and P-wave scattering volumes, in units 
of and m~'^ respectively, are compared with the data of Refs. |42| and |46|. 



length /volume 



Al 



Bl 



CI 



Ref. HI 



Ref. 



ai 
as 
an 
ai3 
asi 
ass 



0.168 
-0.080 
-0.078 
-0.026 
-0.035 

0.201 



0.152 
-0.095 
-0.075 
-0.026 
-0.036 

0.188 



0.142 
-0.105 
-0.070 
-0.025 
-0.035 

0.171 



0.175 
-0.087 
-0.068 
-0.022 
-0.039 

0.209 



0.173±0.003 
-0.101±0.004 
-0.081±0.002 
-0.030±0.002 
-0.045±0.002 

0.214±0.002 
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